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ABSTRACT 

We report the detection of a stellar density cusp and a velocity dispersion increase in the center 
of the globular cluster M54, located at the center of the Sagittarius dwarf galaxy (Sgr). The central 
line of sight velocity dispersion is 20.2 ± 0.7kms~^, decreasing to 16.4 ± 0.4kms~^ at 2"5 (0.3 pc). 
Modeling the kinematics and surface density profiles as the sum of a King model and a point-mass 
yields a black hole (BH) mass of ~ 9400 Mq. However, the observations can alternatively be explained 
if the cusp stars possess moderate radial anisotropy. A Jeans analysis of the Sgr nucleus reveals a 
strong tangential anisotropy, probably a relic from the formation of the system. 
Subject headings: black hole physics — globular clusters: individual (M54) — stellar dynamics 
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1. INTRODUCTION 

The relationships between the mass of central Super- 
massive Black Holes (SMBHs) and the physical prop- 
erties of their host galaxies suggest a close connec- 
tion be tween SMB H and galaxy formation an d evo- 
lution (fKormendy " & Richstonc 1995; Magorria n et al.l 
1998"; 'Gebhardt ct al. 2000; Fcrrar ese fc Merrittl 120001 : 
Graha m et al. 2001). Determining the range of valid- 
ity of these scaling relations and where they break down 
will likely lead to a deeper understanding of the link 
between BHs and their hosts. While numerous ob- 
servational studies have demonstrated the presence of 
SMBHs in massive galaxies, there is only tantalizing ev- 
idence that the correlations also hold in low-mass sys- 
tems, possibly down to globular cluster (GC) scales, 
where BHs with masses of a few 10^ — 10^ M© (com- 
monly named "intermediate-mass BHs" , IMBHs) have 
been sugg ested to reside (see references in iLanzoni et al.l 
120071 and iNovola et all l2008l hereafter N08). For low- 
mass galaxies the properties of the central BH do not 
appear to co rrelate well with host mass Mgai (see, e.g., 
references in lFerrarese et aLll2006l hereafter F06). How- 
ever, F06 have shown that a tight correlation exists 
with the so-called Central Massive Object (CMO), the 
latter being either a BH or a compact stellar nucleus: 

McMO -0.3%Afga/. 

Within this c ontext the system comprising Sgr 
(jlbata et al.lll994D and the GC M54 presents an excel- 
lent prototype to study these issues on different mass 
scales. Sgr is a disrupting dwarf spheroidal, with a con- 
stant surface brightness core of about 1.7 kpc. At its 
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center lies (i) a stellar nucleus (hereafter Sgr,N) made of 
typical Sgr metal-rich stars with a scale-length Vc = 0.4 
pc, and {ii) M54, a v ery bright metal-poor GC (re — 0.9 
pc, M = 2 X 10^ M(T,;'Prvor & Meylan' 1993). As demon- 
strated in Bellazzini et al. (2008, hereafter BOS), these 
two sub-systems are coincident in position and veloc- 
ity, but exhibit different velocity dispersion (cr) profiles, 
thus suggesting that M54 was born in the halo of Sgr 
and was driven to coincide with Sgr,N by orbital decay 
via dynamical friction. Sgr is an excellent candidate for 
a system that would appear to have a compact stellar 
nucleus if it were located at the distance of the Virgo 
cluster. Indeed, the CMO of Sgr appears to be a com- 
bination of the central cluster M54 and Sgr,N, and the 
FOG relation suggests a mass for Sgr of 1 0^ Mp,, in good 
agree ment with pre-disruption estimates (jlbata fc Lewis! 
Il998f ). Here we push the investigation of the validity of 
the FOG relation to even smaller scales, addressing the 
possibihty that M54 also hosts a CMO (an IMBH) in its 
center. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Density Center and Profile 

To measure accurately the density profile and center of 
M54, we re-analysed the high-resolution dataset of B08, 
derived from Advanced Camera for Surveys (ACS) im- 
ages taken in the FGOGM^ and F814W filters, with indi- 
vidual exposures of tcxp = 30 s. To improve the photo- 
metric and astromctric accuracy, the innermost 5" were 
re-measured using ROMA FOT, a package specifica lly de- 
veloped for crowded fields (iBuonanno fc Iannicolal[r 9891. 
The catalog obtained by [Monaco et al.l (|2002l hereafter 
M02) was used outside of the ACS field, while stars in 
common between the catalogs were used to calibrate the 
ACS photometry and astrometry (0''3 uncertainty in ab- 
solute position). 

The densi t y cen ter was measured using the 
iFerraro et al.l ()2003[ ) procedure: from a first guess 
for the cluster center, we computed the mean of the 
positions of all stars brighter than a given limiting 
magnitude and lying within a circle of ~ 8 — 10" 
radius; this value was then used as the new centre for 
the search-circle, and the procedure was iterated until 
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Fig. 1. — (a) ACS color-magnitude diagram (CMD). Sources are classified into likely M54 (heavy black), Sgr (red), and foreground (light 
black), according to CMD location, (b) Corresponding star counts profiles for M54 (black) and Sgr (red), with respect to the M54 density 
center and as a function of (projected) radius R. Towards the centre of M54 we detect a sharp rise in the stellar surface density, at odds 
with a flat distributio n expected from a King model. The green line shows a Hernquist model fit to the Sgr profile (with background given 
by the King model of IMaiewski et aPfSOOSi . also displayed in Fig. 2 of BOS). The background-subtracted profile M54 is presented in (c), 
together with a composite model comprising a King model and a central IMBH. 

convergence. In order to limit contamination we adopted 
the color-magnitude selection criteria shown in Fig. la. 
To avoid biases and spurious effects, we considered two 
samples with two different limiting magnitudes {V < 19 
and V < 19.5). The values thus obtained agree within 
0.2" and we adopted their mean {a = 18'^55'"3.345", 
6 ~ — 30°28'47"1) as the best estimate (cross-hairs 
in Fig. 2 mark this position). Note that this new 
determination is su bstantially di fferent (Aa ^ 0'.'7 , 
A6 5") from the lHarrisI (|1996[ ) center. 

The projected density profile of M54 was measured by 
counting stars down to \^ = 19.5 in 30 concentric annuli 
around this new center. Each annulus was divided into 
an adequate number of sub-sectors; the resulting density 
is the average of the corresponding sub-sector densities, 
and the quoted uncertainty is their rms dispersion. The 
star-counts profile is well matched by the surface bright- 
ness profile (which is unaffected by incompleteness); we 
therefore decided not to apply an incompleteness correc- 
tion, avoiding an additional source of uncertainty. 

The observed density profile of M54 is shown in Fig. lb 
(black dots), where we have also plotted as triangles the 
values obtained by splitting the first bin (r < 0"7) in two 
portions (note that only 6 stars are found at r < 0"35). 
If the innermost (r ^ 1") points are excluded, the density 
profile is well-fit by an isotropic, single-mass King model 
with core radius Tc ~ 6" and concentration c ~ 1.8. How- 
ever, in the inner ~ 1" there is an indication of a power- 
law deviation from a flat core behavior: the profile rises 
as E oc with a ~ 0.3. These values are consistent 

wit h those expected in the presence of a central IMBH 
(see iBaumgardt et al.ir2005l:lMiocchill2007D . Indeed, the 
observed profll e is fit well by a King-|-BH model built 
as described in iMiocchU ()2007[ ). with BH to total cluster 
mass ratio Mbh/M = 0.0047 (i.e., Mbh ^ 94OOM0). 

2.2. Spectroscopy 
The fiber-fed FLAMES spectrograph at the Very Large 



Telescope (VLT) was used to follow up these findings. 
Two observing modes were employed: the "MEDUSA" 
and "ARGUS" configurations, with the high-resolution 



Fig. 2. — (a) Center of the ACS I-band image. The density center 
is marked with cross-hairs, while the pixel with highest a is marked 
with a red circle, (b) Map of a over the same region, derived from 
the 2 deep ARGUS fields. The circles and "star" markers indicate 
the positions, respectively, of M54 and Sgr stars that are brighter 
than 1=16 (triangles mark foreground sources). Bright stars tend 
to lower the measured a in the pixels whose light they dominate. 
Although the density center (0,0) has a high tr, a large region with 
coherently high dispersion is detected to the SW. 

HR21 mode, covering the Ca II triplet. With MEDUSA, 
130 fibers can be allocated over a 25' diameter field, while 
"ARGUS" covers a rectangular 22 x 14 pixel field (ll'(4 x 
7"7). Some observations from July and August 2005 were 
downloaded from the ESO archive, though most of the 
data were obtained on June 2-5 2008. All the data were 
reduced in a homogenous manner using the ESO pipeline 
"Gasgano" software. A total of 944 isolated stars were 
observed with MEDUSA at radii between 0;5 and 12;5, 
while the entire inner region of the system out to a radius 
of approximately 0^6 was covered with 51 overlapping 
ARGUS tiles, each of toxp = 300 s. A further 2 deep 
ARGUS fields (^cxp = 2535 s) were secured on the very 
center of M54/Sgr under excellent seeing conditions (~ 
0"6). The astrometry of the ARGUS data was refined 
by cross-correlating against the I-band ACS image. 

All the individual MEDUSA and ARGUS fibers were 
sky-subtracted and cross-correla ted against an artificial 
template as in Ibata et al.l ()2005D . Velocity uncertainties, 
estimated from the r.m.s. scatter in Gaussian fits to each 
Ca II line independently, are typically below 1 kms~^ for 
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Fig. 3. — The cr profile measured from ARGUS spectra co- 
added in annuli (black triangles and squares) rises rapidly towards 
the cluster center. Beyond ~ 20", the MEDUSA and ARGUS 
spectra can be used to measure the velocities of individual stars, 
uncontaminated by neighbors. The blue and red points show the 
corresponding a profiles for the M54 and Sgr samples, respectively. 
For comparison, we show the earlier BOS results. The line shows 
the kinematic behavior of the King+BH model shown previously 
in Fig. Ic. 

bright stars. 

Towards the center of M54 our limited spatial reso- 
lution causes the superposition of many stars under a 
single fiber. The dispersion in velocities of these stars 
results in a broadening of the observed spectral lines, 
which can be used to measure the stellar kinematics. 
W e used the Penalized Pixel-F itting (pPXF) method 
of ICappellari &: EmsellernI (|2004f ) to perform the line of 
sight velocity distribution deconvolution. Bright isolated 
stars, observed with the same instrument, were selected 
to serve as spectral templates. 

2.3. Velocity dispersion profile 

Fig. 2b shows the a map derived from the two deep AR- 
GUS fields using the pPXF algorithm. The map must be 
interpreted with caution, since we arc in a regime where 
some pixels are dominated by one or a few bright stars. 
Clearly, where a single star dominates, a will tend to be 
abnormally low; furthermore, spuriously large a values 
can be measured in pixels containing only a few stars 
with velocities on opposite tails of the velocity distribu- 
tion. However, the velocity distribution will tend towards 
more reliable values when spectra are co-added over sev- 
eral spatial pixels. 

As is apparent from Fig. 2b, there is a noticeable cor- 
relation between bright stars and those pixels with low a 
values. However, one also perceives a radial decrease in 
a. Surprisingly, the highest a peak (at [—4", — 1"5]) does 
not lie at the density center of the system. Although this 
peak may be an artifact of a chance velocity configura- 
tion of the 4 bright stars that surround it, it nevertheless 
lies in the middle of a ~ 3"- wide high a region, so the 
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Fig. 4. — (a) An MCMC scheme was used to find the most likely 
profiles for the radial (red) and tangential (green) components of 
cr (as a function of radius r) of the Sgr,N population that give rise 
to a projected cr that is fiat as a function of projected radius R. 
(b) The components of cr required of the M54 population if no BH 
is present. (Data points are reproduced from Fig. 3). 

reality of the peak cannot easily be dismissed. Inspection 
of the ACS image (Fig. 2a) reveals that this "kinematic 
center" lies outside any obvious density peak. We ver- 
ified that this does not correspond to the center of the 
metal- rich Sgr stars. 

Bearing this ambiguity in mind, we present the a pro- 
file over the entire radial range of M54 in Fig. 3, with 
panels (a) and (b) obtained with respect to the den- 
sity and kinematic centers, respectively. Circles and 
lozenges mark the a values derived from isolated stars 
observed in single fibers; these values and their uncertain- 
ties were calculated using a maximum likelihood method 
(see BOS). Squares and triangles are derived from AR- 
GUS spectra co-added in annuli, and measured with the 
pPXF software; the corresponding un certainties were es- 
timated by "bootstrap-resampling" (jPress et al.l Il992f) 
the spectra, and re-measuring a in 100 random simu- 
lations at each radial bin. ARGUS pixels within 0"7 of a 
star of magnitude 7 = 15 were rejected. The line shows 
the prediction of the King-|-BH model presented in Fig. 1, 
and gives a close representation of the measured a profile. 

The a profile of the Sgr population is also shown in 
Fig. 3. The isothermal behavior of this population seen 
by BOS at larger radii clearly continues inwards to i? 
20". 

2.4. Anisotropy of Sgr,N 

The data presented above allow us to address a conun- 
drum raised in BOS: how is it possible for M54 (with a 
rising a profile towards the centre) to co-exist in equilib- 
rium with the Sgr,N population, that exhibits an isother- 
mal projected profile, almost constant at cr = 9.2 kms^^? 
To answer this question, we turn t o the Jeans equation 
for a spherically-symmetric system (jBinnev fc Tremaind 
[19871 Eqn 4-55): 



GM{r) 



dlnp 



d In r d In r 



V ad: 



(1) 



where M{r) is the cumulative mass inside radius r, p is 
the density, and cr^ and ag are the radial and tangential 
velocity dispersions, respectively. Interestingly, given the 
Sgr density profile shown in Fig. lb, if we assume that 
the Sgr,N population is isotropic, the total mass profile 
resulting from the Jeans equation in the radial range 1- 
30" turns out to be substantially lower (by a factor up 
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to ~ 10) than that of the King model (Figs. Ic and 3a) 
that best fits the M54 population alone. Since this is 
physically absurd, we conclude that the Sgr,N population 
cannot have isotropic orbits. 

We can now turn this argument around, to ask what 
level of anisotropy is consistent with the M54 mass 
model, since the cluster likely represents the domi- 
nant component in the center of the M54/Sgr system 
a.t R ^ 100" (Sgr stars and accompanying dark matter 
should be relatively unimportant in this region). Hence, 
by assuming the M54 King model for M(r), and the 
Hernquist model fit to the Sgr population for p, ar and 
as remain as the only unknowns in the Jeans equation. 
However, ar and ag are constrained by the observed (pro- 
jected) dispersion a. A trial profile for ar{r), defined 
by a bi-cubic spline at 5 logarithmically-spaced points 
(red line in Fig. 4), was improved upon iteratively us- 
ing a Markov-Chain Monte-Carl o (MCMC) sch eme with 
"parallel-tempering" chains (see lGregorvll2005l ): the set 
target was to produce a flat projected a{R) = 9.2 kms~^. 

Fig. 4a shows that the nested system can be con- 
structed, but it requires a high degree of anisotropy for 
Sgr,N, with orbits being exceedingly tangential between 
20-200" (interior to 20" the fit is based only an extrap- 
olation of the a profile of Sgr,N). If M54 decayed into 
the center of Sgr, as suggested by BOS, its orbital energy 
must have been transferred to stars previously present 
in that region. It is tempting to assume that that this 
would manifest itself as a high ag, as inferred. 

3. DISCUSSION AND CONCLUSIONS 

The star-counts profile of M54 (Fig. Ic) shows a sharp 
rise in the inner 1", a behavior that is also mirrored 
in the a profile (Fig. 3a). Although the rise up to 
a ~ 18kms~^ a.t R — 2" is a robust result, the kine- 
matic data are not as clear-cut interior to this radius. 
It is possible that the stochastic distribution of stars in 
both projected position and line of sight velocity could 
give rise to spuriously high a values. Indeed, we judge 
the "kinematic peak" to be such an artifact, since it 
appears to have no density counterpart. However, the 
central rise in a around the density center is in agree- 
ment with expectations from the stellar density cusp, 
and since this region is relatively free of very bright stars 
that could bias the a measurements, we judge that the 
measured rise in a is likely real. Further support for 
the reality of the M54 a peak comes from a comparison 
to oj Centauri. Recently, N08 analyzed a data-set com- 
prising ACS photometry and Gemini integral-field spec- 
troscopy, finding evidence for both a power-law increase 
in the central stellar profile, as well as a marked rise to 
a = 23kms~^ in the center, which they interpret as pro- 
viding evidenc e for a 4 x 10^ Mq BH. A reanalysis by 
iMiocchil (I2OO9I) finds Mbh/M ~ 0.0057 ± 0.001, which 
(depending on the adopted cluster mass) brings down 
Mb h to 1/7 to 1/2 of the N0 8 value, consistent with 
the Ivan der Marel fc Andersonl (2009) 18000 M© upper 
limit. The similarity to the results presented above is 
striking, all the more as it suspected that lo Cen is the 
stripped core of a galaxy that was once similar to Sgr 
(jBekki fc Freemani r2003. B08). While it could be argued 
that the central dispersion increase in one of these sys- 
tems could be due to chance stellar alignments, it would 
require a conspiracy for both to be affected in the same 



way (especially since the discreteness of the underlying 
stellar distribution should in general lower a) . 

What is the nature of this central structure? Is it 
due to a cold dark matter (CDM) cusp, a population 
of massive stellar remnants, peculiar stellar kinematics, 
or indeed an IMBH? For any reas onable total mass and 
concentration parameters of Sgr, a iNavarro et al.l ()1997f ) 
CDM cusp would possess only ~ 10^ M© within the inner 
10 ", and so cannot ac count for the a rise. As discussed 
in lNovola et al.l ()2008D . neither can massive stellar rem- 
nants be accommodated, because of the implausibly high 
central concentration required if the cluster has not un- 
dergone core-collapse (like M54 and uj Cen). 

In order to explore the alternative that orbital 
anisotropy is the cause of the observed high central a, 
we performed a Jeans analysis similar to that described 
above. We assumed the stellar density profile shown in 
Fig. Ic and the corresponding mass profile (leaving to- 
tal mass as a free parameter). As before, we iterate on 
an initial guess for ar{r) (taking as a prior that M54 
should be isotropic beyond 10"), to find the projected 
a(R) that is consistent with the kinematic observations 
of M54. The results are shown in Fig. 4b: significant, but 
by no means extreme, radial anisotropy (arjag = 1.25 
at 1") is required in the central 2" to reproduce the ris- 
ing a. Rejecting a single strongly outlying data-point 
(at 24"), the reduced-x^ of this model (x^ = 2.7) is 
marginally worse than that of the King-|-BH model of 
Fig. 3a (x^ = 2.5). While very short relaxation times 
in the centre of a cluster are commonly believed to lead 
to isotropic stellar orbits, this might not hold if there is 
a strong density gradient (i.e. if there is a cusp). The 
IMBH solution appears the only viable alternative, how- 
ever, if the orbits must be isotropic. 

As discussed in Macc arone fc Serv illat (2008'), stellar 
winds should provide gas to the central BH and make it 
emit as a radio or X-ray source. Due to the possible con- 
tribution of Sgr gas, such emission might be more signifi- 
cant in the case of the M54/Sgr system, compared to nor- 
mal GCs. However, the expected radio and X-ray fluxes 
are still highly uncertain. Whi le no suitable radio obs er- 
vations are available for M54, iRamsav fc Wul (|2006al lbl) 
have studied its X-ray emission using the Chandra Ob- 
servatory. We retrieved their data from the archive and 
compared the position of the X-ray bright sources to the 
ACS image. Interestingly, with the default astrometric 
solution provided by the Chandra pipeline, their object 
#2 lies within 1" of the density center of M54. Since 
the Chandra absolute astrometric accuracy is 0"6^, while 
the systematics in the absolute position of the ACS are 
~ 0"3, it is possible that this source may be associated 
with the stellar cus p we have identified. Fro m their lumi- 
nosities and colors, iRamsav fc Wul (|2006a[ ) suggest that 
the 7 sources they detected in M54 are likely cataclysmic 
variables or low-mass X-ray binaries. However, their ob- 
ject #2 stands out as being the only irregularly-shaped 
source in this small sample, and with an X-ray luminos- 
ity of Lx = 0.72 X lO^'^ ergs"^ is bright enough to lie 
at the lower limit of the iMaccarone fc Servillatj ()2008f ) 
predictions. 

The data discussed above do not allow us to clearly 
state whether the putative IMBH belongs to M54 or to 

[http : //cxc ■ harvard ■ edu/cal/ASPECT/celmon/ 1 
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Sgr. However, we note that the estimated BH mass is 
consistent with the (extrapolation of the) F06 relation 
only if its host system is M54 since Mbh/Mm54 ~ 5%, 
while a BH 1000 times more massive would be needed 
for a similar mass ratio with respect to Sgr. We there- 
fore may be in the presence of a very interesting system 
composed of a dwarf galaxy hosting a prominent stellar 
nucleus, itself hosting a central IMBH, and both follow- 
ing the F06 relation accurately-verified at higher mass 
regimes. 

If the presence of an IMBH can be confirmed, it would 



have far-reaching implications not only for the forma- 
tion of SMBHs and galaxies, but also for a deeper com- 
prehension of the unclear nucleation process. However, 
as we have shown, the observations can be equally well 
explained by a centrally-concentrated stellar population 
with radially-anisotropic orbits. It will therefore be very 
useful to explore with N-body simulations whether or- 
bital anisotropy can be maintained in such a dense sys- 
tem, and to obtain higher-spatial resolution spectroscopy 
to probe closer to the central mass. 



REFERENCES 



Baumgardt, H., Makino, J., & Hut, P. 2005, ApJ, 620, 238 
Bekki, K., & Freeman, K. C. 2003, MNRAS, 346, Lll 
Bellazzini, M., ct al. 2008, A J, 136, 1147 

Binncy, J., & Trcmaine, S. 1987, Galactic Dynamics. Princeton 

Buonamio, R., & lannicola, G. 1989, PASP, 101, 294 

Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 

Ferrarese, L., et al. 2006, ApJ, 644, L21 

Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Ferraro, F. R., Possenti, A., Sabbi, E., Lagani, P., Rood, R. T., 

D'Amico, N., & Origlia, L. 2003, ApJ, 595, 179 
Gebhardt, K., ct al. 2000, ApJ, 539, L13 

Graham, A. W., Erwin, P., Caon, N., & Trujillo, I. 2001, ApJ, 563, 
Lll 

Gregory, P. C. 2005, Baycsian Logical Data Analysis for the 

Physical Sciences. Cambridge University Press 
Harris, W. E. 1996, AJ, 112, 1487 

Ibata, R., Chapman, S., Ferguson, A. M. N., Lewis, G., Irwin, M., 

& Tanvir, N. 2005, ApJ, 634, 287 
Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, Nature, 370, 194 
Ibata, R. A., & Lewis, G. F. 1998, ApJ, 500, 575 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Lanzoni, B., Dalessandro, E., Ferraro, F. R., Miocchi, P., Valenti, 

E., & Rood, R. T. 2007, ApJ, 668, L139 



Maccarone, T. J., & Servillat, M. 2008, MNRAS, 389, 379 
Magorrian, J., et al. 1998, AJ, 115, 2285 

Majewski, S. R., Skrutskie, M. F., Weinberg, M. D., & Ostheimer, 

J. C. 2003, ApJ, 599, 1082 
van der Marel, R., Anderson, J. 2009, astro-ph/0905.0638vl 
Miocchi, P. 2007, MNRAS, 381, 103 
Miocchi, P. 2009, submitted 

Monaco, L., Ferraro, F. R., Bellazzini, M., & Pancino, E. 2002, 
ApJ, 578, L47 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Noyola, E., Gebhardt, K., & Bergmann, M. 2008, ApJ, 676, 1008 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 

1992, Cambridge University Press 
Pryor, C, & Meylan, G. 1993, Structure and Dynamics of Globular 

Clusters. Proceedings of a Workshop held in Berkeley, 50, 357 
Ramsay, G., & Wu, K. 2006a, A&A, 447, 199 
— . 2006b, A&A, 459, 777 



